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Abstract 

> I 

qs, • We discuss some perturbative techniques suitable for the gauge-invariant treatment of the 
scalar and tensor inhomogeneities of an anisotropic and homogeneous background geome- 
try whose spatial section naturally decomposes into the direct product of two maximally 
q\ \ symmetric Eucledian manifolds, describing a general situation of dimensional decoupling in 
which d external dimensions evolve (in conformal time) with scale factor a(rf) and n internal 
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Ph! dimensions evolve with scale factor b(r]). We analyze the growing mode problem which typi- 
cally arises in contracting backgrounds and we focus our attention on the situation where the 
amplitude of the fluctuations not only depends on the external space-time but also on the 
^ ■ internal spatial coordinates. In order to illustrate the possible relevance of this analysis we 
compute the gravity waves spectrum produced in some highly simplified model of cosmolog- 
ical evolution and we find that the spectral amplitude, whose magnitude can be constrained 
by the usual bounds applied to the stochastic gravity waves backgrounds, depends on the 
curvature scale at which the compactification occurs and also on the typical frequency of the 
internal excitations. 
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1 Introduction 



The assumption of isotropy and homogeneity of the background manifold permits a consis- 
tent theoretical treatment of the space-time evolution of its inhomogeneities which can be 
classified in scalar, vector and tensor modes with respect to the three dimensional spatial 
coordinate transformations on the constant time hypersurface [p], [2|. The different modes 
are decoupled to first order in the amplitude of the fluctuations and this allows the defini- 
tion of perturbed quantities which are invariant under the gauge group of the infinitesimal 
coordinates transformations H |I3|. Within the Bardeen's gauge- invariant approach, the 
amplified primordial spectrum of fluctuations can be reliably computed for a wide class of ho- 
mogeneous and isotropic cosmological models || Q and in particular for the "slow-rolling" 
scenarios leading to the de-Sitter like inflation fTj. 

One of the main motivations in order to relax the assumption of isotropy of the back- 
ground geometry comes from the models of early universe (like superstring theories ||) 
describing the unification of gravity with gauge interactions in a higher dimensional mani- 
fold |J . We will then consider a homogeneous and anisotropic manifold which can be written 
as: 

ds 2 = g^dx^dx" = a 2 (rj)dr] 2 - a 2 (r})j ij dx i dx j - b 2 (r])>y ab dy a dy b (1.1) 

(conventions: fi, u—l,...,D — d + n + 1; i, j=l,...,d; a, b—d+l,...,n; rj is the conformal time 
coordinate related, as usual to the cosmic time t = J a(r))dr) ; j%j(x), jab(y) are the metric 
tensors of two maximally symmetric Euclidean manifolds parameterized, respectively, by the 
"internal" and the "external" coordinates {x 1 } and {y a })- This metric describes the situa- 
tion in which the external dimensions (evolving with scale factor a(r))) and the internal ones 
(contracting with scale factor 6(77)) are dynamically decoupled from each other. In order 
to compare the phenomenological consequences of the models formulated with extra dimen- 
sions it seems crucial to correctly compute the amplified spectrum of inhomogeneities but, 
unfortunately, the treatment of the fluctuations in an anisotropic manifold becomes quite 
cumbersome also because of the natural coupling arising among scalar vector and tensor 
modes to first order in the amplitude of the fluctuations. The original investigations in the 
subject H stressed that the discussion of the fluctuations in the synchronous gauge is com- 



1 



plicated also because of the spurious gauge modes already present JL4| in the isotropic case. 
The problem of the metric fluctuations in an anisotropic background geometry was then 
addressed within the Bardeen formalism with two different and complementary approaches 
nj. It was actually shown that the gauge invariant quantities can be constructed not 



only, separately, in the external and in the internal manifold |]15 |, but also over the whole 
manifold of Bianchi-type I Jl6| (gravitational waves in Bianchi-I universes were also discussed 
/J). Following the first of the two mentioned approaches it is possible to distinguish, 



in 



from the purely mathematical point of view, the scalar vector and tensor modes in each of 
the two manifolds. Even though this classification will be technically very useful it does 
not necessarily coincide with the physical situation (as correctly stressed in since, for 

example, the tensor fluctuations polarized along the internal dimensions will be seen by 
an observer living in the external space as scalar fluctuations. The evolution equations for 
each type of perturbations were then solved well outside the horizon under the assumption 
that the Laplacians belonging to the external and internal manifolds were negligible. The 
anisotropic extension of the scalar Bardeen potentials were shown to grow much faster than 
the tensor and vector gauge-invariant amplitudes in the vicinity of the collapse of the internal 
scale factor. The very fast growth of the scalar modes outside the horizon was also discussed 



in the context of the dilaton-driven solutions in string cosmology [18|, where it was found 



24j that even though the rate of increase of the scalar fluctuations is much faster than in 



the usual inflationary models characterized by a quasi de-Sitter spectrum || a perturbative 
treatment is still plausible, at least in the (3 + l)-dimensional case with static internal dimen- 
sions, by carefully "gauging-down" the scalar growing modes. In spite of these attempts, in 
an anisotropic background the solution of the evolution equations well outside the horizon 
does not suffice, by itself, for the calculation of the spectrum of metric perturbations. In 
order to give a reliable expression for the space-time evolution of the proper amplitude of 
the fluctuations it can be assumed |2]J relying on the particular features of the background 
evolution, that the only effective dependence of the perturbed quantities from the internal 
dimensions comes in through the time evolution of the compactification radii. This approxi- 
mation scheme can be illustrated using the evolution equation of the tensor modes polarized 
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along the external dimensions (which we will derive in Sec. 3): 



hi" + [(d - l)H + nf]h( - Vlhi - a -^\ = (1.2) 

(' = d/dr/; Tt = (In a)', T = (In b)'\ V|, and V| are, respectively, the external and the internal 
Laplace-Beltrami operators). In some particular model it can actually happen that a<b, at 
least for scales which went out of the horizon before the compactification was achieved. In 
this case it is possible to neglect the internal gradients compared to the external ones, and 
all the dependence from the internal dimensions will be given by T which vanishes in the 
case of static internal scale factors. If, on the contrary b<a, the internal Laplacians cannot 
be neglected especially prior to the dimensional decoupling when the internal and external 
scale factors were of the same order. Since the tensor modes only couple to the background 
curvature in order to discuss their evolution we only need to specify the time evolution of 
the scale factors. The scenario we want to examine consists in general of two phases. A 
multidimensional phase where the time evolution of the scale factors can be parameterized 
as 

a{rf) ~ \r]\ a , 6(77) ~ \rjf . (1.3) 

This phase can be genarally followed by a compactification phase which glues together the 
multidimensional epoch and the ordinary, isotropic, FRW evolution. There are different com- 
pactification scenarios corresponding to the parameterization (|1.3|) . In It is actually possible 
either to assume that the compactification occurs at some stage after the initial "big-bang" 
singularity (i.e. for 77 > in (|1.3j)) as, for instance, in |10|, |25[ ) or at some stage before the 
"big-bang" singularity (i.e. 77 < in ( |1.3j )) as seems more likely in the pre-big-bang scenarios 
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lq, |19| . Even though our considerations, at this stage, are purely kinematical, it is anyway 



useful to point out that the two mentioned compactification pictures are dynamically very 
different since in the pre-big-bang case the end of the higher dimensional phase could coincide, 
in principle, with the beginning of the ordinary isotropic evolution while in the first picture 
some other mechanism is required in order to smoothly connect the multidimensional phase 
(trapped among two singularities) to the FRW universe. An interesting issue is then if in 
the context of the string inspired models of cosmological evolution the usual problems of the 
ordinary Kaluza-Klein models (stabilization and isotropization of the internal dimensions, 
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backreaction effects due to particle production |42|, [43]]) can be solved (or at least alleviated) 
by the mechanisms usually proposed in order to regularize the time evolution of the 
curvature invariants and in order to slow down the dilaton growth. Since at the moment this 
solution is still unclear not only in more than four dimensions but even in some simplified 
two-dimensional toy model of cosmological evolution |27j] (where the quantum backreaction 
was taken into account and where the problem of the dilaton seem growth seem still to per- 
sist) we will concentrate our attention on the main kinematical features of the dimensional 
decoupling. Our purpose will only be to stress the regimes of the time evolution of the scale 
factors where, possibly, the internal Laplacians of eq. ( |1.2|) are leading if compared to the 
internal ones. More specifically for t > an eventual accelerated expansion in the external 
space (a > 0, a > 0, i. e. a < — 1 in ( |1.3| )) required in order to solve the kinematical 
difficulties of the Standard Model, together with a simultaneous contracting evolution of the 
internal dimensions (b < 0, (3 < in ( |1.3|) ), would not forbid in the large t limit — > 
if a < — 1) the dominance of the internal Laplacians over the external ones (a < (3 in our 
parametrization). Similar conclusions can be reached in the limit 77 — ^ — if, for t < 0, we 
consider an accelerated contraction in the external space (a < 0, a < 0, — 1 < a < 0) as 



suggested by the Einstein frame picture of the string cosmological models [p3|, |25| in order to 
solve the flatness and the horizon problems. These qualitative considerations suggest that 
if a given external Fourier mode k went out of the horizon (kr] ~ 1) during an early phase 
where a> b, the contribution of the internal Laplacians have to be seriously considered. At 
the same time a complete solution of the evolution equations of the scalar and tensor fluctu- 
ations depending on the internal and external spatial coordinates was never studied not even 
in some oversimplified model of background evolution. A reliable computation of the power 
spectrum for scalar and tensor inhomogeneities in higher dimensional theories is beyond the 
scope of the present investigation, nonetheless we would like to analyze the evolution and 
the amplification of the metric fluctuations in some specific toy model with extra dimensions 
but without assumptions for what concerns the evolution equations of the fluctuations. We 
would like also to avoid any kind of "slow-rolling" hypothesis in the solution of the back- 
ground equations which could confuse the analysis of the perturbations. For this reason we 
shall mainly discuss two classes of exact solutions of the multidimensional Einstein equations 
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: the Kaluza- Klein "vacuum" solutions which can represent a good description in the vicinity 
of the collapse of the internal dimensions and the multidimensional anisotropic universe 
filled with scalar field matter. We find quite useful to work from the very beginning with 
the following scalar-tensor action: 

S = S 9 + S m = ~ J d D x^R + J d D xV=g{~g af3 d a( pdp<p - V(<p)} (1.4) 



(where Id = ySirGo/S', if V = and tp = the action can describe a vacuum Kaluza-Klein 
phase and if <p' 2 » a 2 V we recover the tree level string theory effective lagrangian [in 
D = 10 critical dimensions] for the massless modes of the theory, written in the Einstein 
frame and in the absence of antisymmetric tensor field) . 

The plan of the paper is the following. In Sec. 2 we will review the higher dimensional 
background equations of motion and the particular classes of exact solutions which will be 
used in the following sections as theoretical laboratory for the analysis of the fluctuations. 
In Sec. 3 the Bardeen approach for the scalar and tensor perturbations will be discussed. 
Particular attention, in the case of the scalar fluctuations, will be paid to the possible gauge 
choices which completely fix the coordinate frame and to the diagonalization of the system 
of perturbed equations. In Sec. 4 we will focus our study on the evolution of tensor pertur- 
bations and we will compute the normalized spectral amplitude for two simplified models of 
dimensional decoupling. In Sec. 5 we will move to the analysis of the scalar inhomogeneities 
and we will approach the growing mode problem within the formalism discussed in the 
previous Sections. Sec. 6 contains few concluding remarks. 

2 Background models 

The variation of the action ( |1.4|) with respect to g^ v and to <p provide the equations of motion 
for the background fields: 

R u , - = 3^t; (2.1) 

dV 

9 aP V a V^ + — = (2.2) 
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(where = d il ipd v ip — gfa,{\g a ^ d a (pdp<p — V((p))). If we restrict our attention to the case in 
which the scalar field is homogeneous (ip = ip(n)) the evolution of the geometry is completely 
determined by the time evolution of the two scale factors a(i]) and b(n). 
Using the line element ( P-.1Q , eqs. (|2.1|) - ([2.2|) become: 



d(d- 1 



n(n — r 



H + — -T 2 + ndHT 



d(d-l) n(n-l) a 2 

^ l^a H ^ IT^b 



J2 



b 2 



3l 2 D (^- + a 2 V), 



(00) 



(d - l)H> + ^'^'V + nT> + ^ 1} - 2 
, {d-2){d-l) ir n {n-\) a 2 ^ 

H " A-a H " TT^ 



F 2 + n{d-2)HJ : 

/2 



b 2 



— w 2 + -^- r 



3l 2 D (a 2 V-r-) 



{n - + dH' + 



,d(d-l) ( w - 2)^-1) ^ 
H - H ~yz b 

dV 

dip 



1) ^ 



^ 2 + (rf-l)(n-l)W 
) (aa) 



3l 2 D (a 2 V - r- 



(2.3) 



2 2 

+ [(d - 1)H + nf)<p' + = ( r - ) 

(where K, a and /C& are respectively the curvatures constants of the external and internal 
maximally symmetric spaces). In this paper we will generally work in the case of spatially 
flat internal and external manifold with topology Md+i <8> T n (where d+ 1 is the conventional 
(c?+l)-dimensional flat universe and T n is an n-dimensional torus). Summing and subtracting 
the previous equations we get two useful relations: 

r = -T\nT+ (d- 1)H] + 



-U\nT + {d - l)H] 



6l 2 D a 2 V 
(n + d-1) 
6l 2 D a 2 V 



(2.4) 



(n + d-1) 

If V — and ip = the solutions of the system (|2.3|) define the Kaluza-Klein vacuum 
A class of exact vacuum solutions can be obtained using the power law ansatz of eq. (|L 
in the equations of motion (|2.3p 



a 







d±Jd 2 + d(n + d)(n-l) 



d(n + d - 1) T Jd 2 + d(n + d)(n - 1) 



nd =F dJ d 2 + d(n + d) (n — 1) 



nd(n + d — 1) =F nJd 2 + d(n + d){n — I) 



(2.5) 



(the exponents of the scale factors in cosmic time (a = a/(a + 1), (3 = j3/(a + 1)) are related 
by the Kasner sum rules <i<S + n/? = 1 and da 2 + nf3 2 = 1). The twofold ambiguity in the sign 
of the exponents shows that there are two independent solutions for each number of internal 
and external dimensions. A particularly simple case which will be used in our analysis is the 
solution with n = 1 and d = 3. In this case the two solutions are a — 0, /3 — 1 and a = 1, 
(3 = — 1. Since we want to analyze mainly the contribution of the internal dimensions to 
the evolution of the fluctuations we choose the solution i which the external dimensions are 
static (a = 0, j3 = 1) so that all the contribution to the amplification of the fluctuations will 
come, effectively, from the internal space. The contracting branch (/3 = +1, i. e. 6 = and 
6 < since t ~ if) can be matched with a radiation phase 

a{V) = 1, Kv) = ( ), V < ~Vc 

Vc 

a( V ) = (^^), 6(77) = 1, 7]>-r, c . (2.6) 

Vc 

This toy model is not realistic and somehow artificial since the radiation is not dynamically 
generated but only assumed. In a more refined treatment the back-reaction effects should be 
correctly taken into account j|2], |45| since the scalar and tensor inhomogeneities amplified 
during the classical evolution can eventually modify the background dynamics leading to 
an effective damping of the anisotropy of the background metric (as usually happens in the 
(3 + l)-dimensional anisotropic models of Bianchi type-I |44[]). Nonetheless (|2.6|) shares some 
essential features of a realistic scenario of dimensional reduction in which the (3 + 1) external 
dimensions decouple from the fifth one down to a compactification scale H c ~ l/rj c and our 
purpose will be to connect the amplitude of the scalar and tensor fluctuations not only with 
the curvature scale but also with the typical frequency of the internal oscillations (which can 
be also constrained, with different arguments, from the present value of the fine structure 
constant flCf). 

In the case of negligible potential, using the previous power-law ansatz for the scale 
factors from eq. (|1.2|) and a logarithmic ansatz for the scalar field 

¥ ~ ~~ t= h 1 \v\ (2-7) 
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another class of solutions of eq. (|2.3| ) is given by: 



=F(1 - d — n) ± n =F (d — 1) — y/d + n 
a = 



(Vd + n + l)(d + 



n 



±(1 — d — n) ± n =p — 1) — \J d + n 
(VdTn+ l)(d + n - 1) 



_ / 2 ±ti T (j - 1) - Vd + ^ - 1 

By choosing everywhere the upper sign and in the case of critical dimensions (D = 10 with 
d = 3 and n = 6) the solution fl2.8Q is a particular example of the general exact dilaton-driven 
solutions originally derived |24| in the String frame usually related to the Einstein frame by a 
conformal rescaling of the metric tensor (g^l rtn9 = g^ nstem exp [(2ip — 2(pi)/ (d + n — 1)]; the 
dilaton is only redefined according to tp Einstein = \j^/{d + n — i)<f string Wd = 1 so that ip is 
dimension- less] ) . We will discuss, as a particular toy model, a 10-dimensional dilaton-driven 
solution continuously matched with the radiation phase: 



< -7 



c 



a iv) = ( Kv) = 1 3 V 2 = const , r/ > -r/ c . (2.9) 

% 

In this model for r/ < — r/ c we have an accelerated contraction of the external dimensions 
supplemented by the a decelerated contraction of the internal ones. We also notice that for 
7] < —r] c the scale factor are related by a duality relation (a = 1/b) which more generally 
holds in the String frame (a = 1/b) provided we choose the upper sign in ( |2.8| ). We point 
out that this model is not realistic for the same reasons mentioned in the case of ( |2.6| ) and 
also because it was shown that in order to have a graceful exit from the dilaton-driven epoch 
it is crucial to include in the picture a stringy phase during which the background dynamics 
is driven by the higher order in the string tension expansion [25|| . In both the examples ( |2.6| ) 
and (|2.9|) the internal scale factors are static during the radiation dominated era, while a 
time dependence could be, in principle, also included in the internal scale factors during the 
radiation and matter dominated epochs. The time dependence in the internal scale factors 
would be anyway strongly constrained by nucleosynthesis which would require pi , |52]] during 



the radiation dominated epoch b ns /bo < 1 + e (b is the actual value of b and |e| < 10 



-2> 
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while in the matter dominated epoch the constraint would be instead b/b = jF/a < 10~ 9 H 
(H = 1.1 x 10~ 28 /iioo cm" 1 ). Moreover the time variation of Go is also constrained P3"| 
since G nijiucl) j 'G d(vo) = 1 + e (|e| < 3 x 10 _1 ) in the radiation dominated epoch, and 
\Gd/Ge>\ < 10~ 1 H during the matter dominated epoch. In our naive models 6=1 (for 
rj > —rjc) the previous constraints are automatically satisfied. 



3 Scalar and tensor fluctuations 

The scalar and tensor fluctuations of the geometry ( |1 . 1| ) can be discussed within a gener- 
alization |Tj| of the gauge- invariant formalism || [Tj|]. The infinitesimal coordinate trans- 
formations preserving the scalar nature of the fluctuations with respect to each maximally 
symmetric space are: 

x l —> x 1 = x % + e l 

r] — > fj = T] + e° , (3.1) 

(where e° = e (x,y,r]), ( a = d a ((x,y,r]), e % = d l e(x ) y,ri)). The perturbed scalar metric can 
be written in terms of 8 linearly independent scalar quantities 

/ 2a 2 4> -a 2 Bi -abC a \ 

SgijD = - a 2 Bi 2a 2 ip5 ij - 2a 2 E ij -abD ia (3.2) 

\-abC a -abD ia 2b 2 £5 ab -2b 2 G ab J 

(conventions :Bi = B^, E^ = E\ij, C a = C\ a , G ab = G\ ab , D ia = D\ ia ; the bar denote the 
covariant derivative with respect to one of the two internal spatial metrics depending on the 
index, and it coincides with the ordinary partial derivative if K a — K\, — 0). The fluctuations 
in the scalar field will be 

<f(v,x,y) -> (p(ri) +x(v,x,y) • ( 3 - 3 ) 



Under an infinitesimal coordinate transformation ( |3.1D the perturbed scalar quantities change 
as follows: 



-He 



0' 
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tfj _> ^ = ^p + ne° 

C - £ = Z + Fe° 
E -> E = E — e 
B -> £? = £ + e ° - e' 

c - c = c + V--c / 

o a 

D - £ = D-*C-^ 
a o 

G -> G = G-C 

X -> X = X-<p'e° • (3-4) 

A possible set of linearly independent gauge-invariant quantities is then: 

$ = +![(£-£>]' 

= i/j + H(B-E') 
E = t + HB-E') 

= D - -G - \e 
a b 

X = x + (f >(B-E') . (3.5) 

Notice that \I> and $ coincides, up to a sign, with the Bardeen's potentials while 5, Q and 
appear only in the anisotropic case. In the homogeneous and isotropic case it is always 
possible to choose a particular coordinate system by completely fixing, to first order, the 
arbitrary scalar functions appearing in the transformations (|3.1| ). If the scalar functions are 
completely fixed (like in the case of the conformally newtonian gauge [[| [T^j ) the equations 
of motion of the fluctuations will be second order differential equations, if, on the contrary, 
the infinitesimal scalar functions are not completely fixed (like in the case of the synchronous 
gauge [0, ||, |14|) the evolution equations will be of course linear but of higher order. Since 



we want to make our problem more tractable we completely fix the coordinate system 

D = 0, B = 0, E = 0, (3.6) 
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in eq. ( |3.4|) and, as a consequence, e°, e and ( are determined from the same equation. In this 
gauge the longitudinal fluctuations (0, £, if)) coincide with the corresponding gauge-invariant 
quantities defined in ( |3.5| ) and in this sense it represents a generalization to the anisotropic 
case of the conformally newtonian gauge. By perturbing to first order the Einstein equation 
( pTTD and the scalar field equation ( [2.2|) we obtain: 

5R U » - \^{9 aP 5R aP + 8g^R a8 ) = 3l 2 D 5T» (3.7) 
8J*{d a dptp - Yip.*) + g aP {d a d pX + SV^d a <f + rzpd a x) - = (3.8) 

(ST^g and 8R^ V are the perturbed affine connections and the perturbed Ricci tensors, the 
indices are raised using always the background metric g^ v ). By using the background field 
equations (|2.1|)-(|2~2j) we can write down explicitly the evolution equations of the fluctuations 
given the perturbed form of the metric ( |3.6| ) in the gauge ( |3.4| ). The (i ^ j) component of 
eq. (|3.7p implies 

<j) = (d - 2)i; + n£ - V 2 V G (3.9) 

which allows to eliminate <f> from all the evolution equations of the perturbations. From the 
(Oi), (Oa) and (aj) components of ( |3.7|) we get, respectively: 

(d - 1)1// + (d- 2)ip[(d - l)H + JF] + n? + n£[(n + 1)F + (d - 2)H] 
-{Vp)' - [(d - 2)H + {n+ l)F]Vp + = 31 2 dV ' X (Oi), (3.10) 

- + ( n ~ + + ^ d " V n + ( d - 2 )( n - 

+ « - Vp){dH + (n- = 3l 2 D <f'x (0a), (3.11) 

- ^[(d - 2)H + (n + l)f]C - + if, - £ + V?G = (aj), (3.12) 

which are not equations of motion but only constraints connecting the fluctuations to their 
first time derivative. The repeated use of eq. ( |3.9| ) together with the equations of motion 
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of the background fields allows finally to write down the (00), (i = j), (a — b), (a / b) 
components of the evolution equations 

- \iypT + ~(V§G)'[(n - 2)T + (d + l)H] - -V^G - —VffiG 
+ nV 2 £ - £'[ndH + n(n - 1)7] + (d - l)V 2 ip - ip'[d(d - 1)H + ndT\ 

+ Yb {v2 y c) '~k [dn + {n ~ m ^ c 

= 3l 2 D [<p'x' + 2a 2 V((d-2)^ + nZ-Vp + a 2 V'x)] (00), (3.13) 

- \(y 2 G)" - l -{Vp)'[{M - 5)H + (3n + 2)T] - \v 2 Vp - ~V*Vp 
+ n£" + g[2n(d - 1)H + n(2n + l)F - nH\ 



+ (d - l)ip" + ip'[(d - l)(2d - 3)H + n(2d - 2>)F\ 

2b K 'y ' 2b" ' y ~ 



+ ^(Vp)'-^(Vp)[(d-2)H + (n + m 



= 3l 2 D [<p'x'-a 2 V'x-2a 2 V((d-2)i> + nt-Vp)} = (3.14) 

- \{Vp)'[{M - 1)H + (3n - 2)^ - \v 2 Vp - 

- V 2 ip + #" + V'[2d(d - l)H + 2(d - l)(n - 1)^] 

+ V 2 £ + (n - + f'[(d(2n - 1) - (n - 1))W + 2n(n - 1)^] 
+ |(V|C)[^+(n-l)^] + |(V 2 C)' 

= 3/ 2 ,[^V-aV X -2aV(( C /-2)^ + <-V 2 G'] (a = 6), (3.15) 

a 2 

G" +[(d - 1)H + nT\G' - V?G - — V§G 

- |C-|[dW+(n-l)J=]C + 2(V»-0^ = (a ^6), (3.16) 

The linear system of differential equations with time dependent coefficients formed by the 
three constraints (|3~TH| ), (gTTTJ) and (^JJ)and by the equations ( gig) , ( gig ), (|XT3| ), (151BQ 
determines the classical space-time evolution of the five fluctuations ?/>, £, G, G and m 
order to simplify the system we also write the perturbed equation of motion for the scalar 
field which can be obtained from the combination of the other equations: 

- V\ + X" + [(d - l)H + nF) X ' - 2(f/[(d - IW + n£') 

- V"a 2 X + 2a 2 V((d-2)ij + nZ-Vp)=0 (3.17) 

12 



Subtracting eq. (|3.14|) from eq. (|3.13| ) and eq. ( |3.15|) from (|3.13|) we get respectively: 



n 



□A + 3[(d - l)H + nf\\' - (2(Vp)' - ^p))[H + —7] 
fi/ 2 

\2 V72 



(3.18) 



(where A = ip + g^f , □ = (d/dr)) 2 - V? - fjV?) and 



n — 1 



} 



+ d{nip + tp' 

7 2 



(3dn — d — n + 1) + SnJ 7 
J" 

3(d - l)H + — (2(d - l)(n - 1) + nd) 



(2(V^0'--(V^)[dtt+(n-l)^ 
6/? 



(3.19) 



- tty^Vx + 2a 2 V((d - + < - Vp)} = 0. 

Combining now the constraints ( |3 . 1 2| ) with eq. ( |3. 16| ) we get an useful expression which 
allows to eliminate C from the other equations 



C 



1 



■\G" + G'((d - l)H + nT) - VlG - -V?G] 



2(H-F) l ~ '"^ -x" 6 2-?/~j • ( 3 - 20 ) 

Using ( |3.20| ) in ( |3.18| ) and ( |3.19|) together with the background equations in the form (|2.4|) it 
is possible to show, by linearly combining the obtained relations, that the evolution equations 
for the longitudinal fluctuations can be written in the gauge ( |3.6|) as 

(3.21) 
(3.22) 



w 



4 

V 


-Vlv- 


a 

— V-v 


7" 

— w - 

Z2 


-V 2 s w- 





where 



v = z 2 X + ziX, w = — 



n(n + d — 1) H JF 



and 



a 2 1)2 (p 



Since 0, ip and £ coincide, in the gauge 
tities listed in eq. ( |3.5| ) we can write that 



d-l , n r — , i 

z 2 = a 2 62 = — K 
a 2 



(3.23) 



(3.24) 



with the corresponding gauge-invariant quan- 



V = z 2 X + ziA, A = ¥ 



d 



W 



Zl 



" { " + (1 - 1 ^ 7 is-^). (3.25) 



h\ 6(d-i; 
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Notice that in the absence of internal dimensions W is zero and V coincides with the scalar 
normal mode of oscillation which diagonalizes the the action ( |1.4j ) perturbed to second order 



in the amplitude of the fluctuations | 12| . By solving the equations for V and W and by 
using their definition in terms of the gauge-invariant fluctuations it is possible to obtain the 
time evolution of all the quantities listed in (|3.5|) (the explicit solution for the longitudinal 
fluctuations will be discussed in Sec. 5). 

The evolution equations for the gauge-invariant tensor modes (only coupled to the scalar 
curvature and not to the sources of the background) can be obtained without any specific 
gauge choice. The form of the perturbed metric will be in this case : 

/0 \ 



1/1V 



—a 2 hij 



(3.26) 

VO -b 2 H ab J 
with V{h^ = h] = and V a H ab = H% = (since H h a and h\ are pure tensor modes in each 
space they are also automatically gauge-invariant with respect to gauge transformations 
preserving the tensor nature of the fluctuations in the external and internal manifold). The 
evolution equations can be easily written by perturbing the Ricci tensor. All the components 
of the Einstein equations (|3.7|) are zero but and (a, b): 

h" + [(d - l)H + nT\ti - V%h - ^V 2 yh = 0, 
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H" + [(d - l)H + nT\E' - V%H - t^V?# = 0, (3.27) 

where h\ = h(r],x, y)el and H% = H(rj,x,y)El ( and E b a are, respectively the external and 
internal polarization tensors). As we already reminded the gravity waves polarized along the 
internal dimensions will not be able to excite a detector of tensor waves but will be seen in 
the (d + 1) -dimensional "external" world as density fluctuations and this means, from the 
mathematical point of view, that h\ are scalar eigenstates of the internal Laplace operator 
(V?) in the same way as H h a are scalar eigenstates of the external Laplace operator (V|) (of 
course h\ and H h a are also solution of the tensor Helmotz equation, respectively, in the exter- 
nal and internal manifold). By defining \i = (l^Alr^a^b^h and M. = (l/24l D )a i i l b% H 
equations (|3.27|) can be easily rewritten as : 



^9 . ,// 9 . , 9 . , ^9 



H" - V?/x - V|// - = 0, M" - V%M - V 2 yM - -^M = (3.28) 
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We notice that in the absence of the internal dimensions pi coincide with the amplitude of 
the tensor normal modes of oscillation of the action (|1.4p fll2"| . The solution of the coupled 
systems of differential equations which describe the evolution of the scalar and tensor inho- 



mogeneities was previously investigated [15, |24| well outside the horizon (where the internal 
and external Laplace operators are subleading) in the case of the Kaluza-Klein "vacuum" 
solutions (|2.5| ) and in the case of the dilaton-driven solutions ( [2.8|) . Provided we neglect the 
internal and external Laplacians it can be actually shown quite easily that the solution to 
eq. Q3gf) is : 



M = a ^bHA 1 + Bl j^), M {n) = a^bHA 2 + B 2 j^), (3.29) 



where Bi are the integration constants. Since from ( |2.5| ) and ( |2.8| ) a ^bz = for 



both the backgrounds in an arbitrary number of dimensions, we will have (from eq. ( |3.29D ) 
that outside the horizon the gravity wave amplitude diverges at most logarithmically : 



/i~ ^J\t,\(A 1 + B x In M), M~ y\rj\ (A 2 + B 2 \n \rj\) . (3.30) 

In the case of the scalar inhomogeneities it can be shown instead (for example from eq. 
( p,18| ) neglecting the Laplacians and for V = 0) that the longitudinal fluctuations diverge 



like a power outside the horizon for both the background solutions ( |2.5|) -( [2T8|) and typically 
we will have: 

A~ a + % (3.31) 
rj 

(a similar behaviour can be deduced also for \1/ and H). If growing solutions are present 
it is in general not possible to keep the amplitude of the fluctuations small all the time 
and at some point the perturbative approach will unavoidably break down leading to the 
so called "growing- mode problem" (which will be addressed in Sec. 5). In order to reliably 
compute the spectrum of the scalar and tensor fluctuations it seems important to consider 
explicitly the contribution of the internal and external Laplacians and since this could be 
quite difficult for a generic multidimensional background we will limit our attentions to the 



models described in eq. ( |2.6|) and ( |2.9| ) which will be analyzed in the following two Sections. 
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4 Graviton production from extra spatial dimensions 



The spectral energy density of the cosmic gravitons produced thanks to the time evolution 
of an homogeneous and isotropic cosmological model is an important source of dynamical 
informations @, E| since the slopes of the spectra versus the frequency offer a snapshot 



of the early history of the hubble parameter p4| , p5| . Moreover the graviton spectra can 
indirectly constrain the homogeneous and isotropic inflationary models [36 , |20| . If we relax 
the assumption of isotropy of the background manifold it is unclear how to perform the 
calculation of the spectral amplitudes which will have to be eventually compared to the 
phenomenological bounds (see p(J and [0 for two reviews concerning the stochastic gravity- 
waves backgrounds and their detectability). Our purpose is to calculate the graviton spectra 
produced in the two oversimplified models of dimensional decoupling presented in Sec. 2 
in order to get the feeling of what could happen in a realistic situations. We will see that 
thanks to the coupling among the scalar and the tensor modes the gravity wave evolution 
equation will get a "massive " contribution which might also be relevant in the case of more 
motivated background geometries. In terms of the eigenstates of the Laplace operators 

V 2 Mk,q) = -k 2 h{(k,q) 

Vlh{(k,q) = -q 2 h{(k,q), (4.1) 

eq. dp|) will be 



^" + [k 2 + q 2 ^- Z \ = . (4.2) 
We consider first of all the model fl2.6|) . For rj < —rj c eq. ( |4.2|) becomes : 

// + [k 2 + JL + ^ = 0, (4.3) 
whereas for rj > — rj c the same equation will be 

fi" + [k 2 + ^( V + 2 Vc ) 2 ]fi = 0. (4.4) 



For r\ < — 7] c an exact solution of eq. (|4.3|) can be written in term of the Hankel functions: 

KkV, q) = ^\[kr]H (2 \kT]), v = i{qr} c ) (4.5) 
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(we have chosen the positive frequency mode which corresponds [for rj — > — oo] to the Bunch- 
Davies vacuum). In the absence of the internal Laplacians ( or, equivalently, if we would 
keep only the lowest mode of the internal excitations, q = 0) instead of ( |4.5| ) we would get 
a completely different solution: 

^kr,) = ^=y/hjH^(kn) (4.6) 

whose limit for small arguments holding when the given mode is well outside the horizon 
(kr) « 1) gives : 

_ 2% 

fj,(krj) — y/rj(l \nkrj) (4.7) 

7T 

which is clearly consistent with the evolution of the gravity waves outside the horizon in an 
arbitrary number of dimensions (and for generic initial conditions) obtained in ( |3.30| ) For 
T) > —7] c the solution of eq. ( |4.4j ) can be written in terms of parabolic cylinder functions 



(defining a new variable z = y2q/r] c (r] + 2rj c ) eq. (|4.3| ) becomes exactly the parabolic cylinder 
equation expressed in its standard form []52], [53]]): 



/j(kr), q) = (— )*(c-E(a, z) + c + E*(a, z)) (4.8) 

Vc 

(a = —k 2 i] c /2q and E(a,z), E*(a,z) are complex conjugated solutions of the parabolic 
cylinder equation |53|). If a > z 2 /4 namely if k 2 > q 2 ((r] + 2t] c )/i] c ) 2 we have that the 
solution ( [4.8|) becomes : 

fi{k V , q) -> —( c _ e - ik ^ 2 ^ + c + e lk(v+2r,c) ). (4.9) 
In the opposite limit (k 2 < q 2 ( v+ ^ Vc ) 2 ) solution (|4.8|) becomes instead 



Kkv, q) ^ J q{r] + 27]c) ( c - e + c+e-* q —*-) (4.10) 

(as can be directly obtained by solving eq. ( |4.3|) for a negligible k 2 ). The last solution 
is identical to the evolution equation of a minimally coupled scalar fields in the radiation 
dominated era with mass m ~ q so that the effect of the internal Laplacians on the evolution 
of an externally polarized gravity wave evolving during the radiation dominated era can be 
described with an effective mass term whose magnitude depends on the magnitude of the 
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excitations belonging to the internal space. In the Schroedinger-like equation (|4.2j )- (|4.3|) the 
mass term modifies the potential barrier whose maximum is now l/4r7 c 2 + q 2 . This effective 
potential barrier leads to wave amplification 0, ^J, or, equivalently, to particle production 
||30|| . Actually the positive frequency modes (for 77 — > —00) in eq. Q4.5|) will be in general 
a linear combination of modes which are of positive or negative frequency with respect to 
the vacuum to the right (rj — > +00). The coefficients of the Bogoliubov transformation 
(c±, |c + | 2 — |c_| 2 = 1) connecting the left and right vacuum and fixed by matching, in 
77 = — rj c , each solution and its first derivative will determine the spectral density of the 
produced gravitons. We can now compute the amplification coefficient c_ in the sudden 
approximation |4l|], namely for (kr] c ) 2 < 1 + (qi] c ) 2 (for (ki] c ) 2 > 1 + {qi] c ) 2 there is no wave 



amplification and the Bogoliubov coefficient c_ is exponentially suppressed ). We will have in 
general a two branches amplification coefficient depending if the mode k is "non-relativistic" 
(k 2 < q 2 ( nd ^) 2 ) or "ultrarelativistic" (k 2 > q 2 ( It ^) 2 ). So matching the solutions (gj ) 
and (|4.9| ) in rj = — r] c we obtain: 

e**c_ * - u^bhyu-i _ l{ ^hy-\ {v + I) _J_ ], (4 .n) 

2^2 2 ; vr 2 ; v 2 ; v 2 ; z/I» J ' v ; 



for A; 2 > g 2 ^M)!, and 



e^c_ ~ ^[^+ (1 - (4.12) 

Z7T 2 



for /c 2 < q 2 q ^ 2 ,?c - ) . The typical amplitude of gravity waves over scales k 1 is 5h(k,q,r]) = 

n+d—1 

l D k d/2 q n/2 \h(k,q,r])\ (where l D = M P 2 ) and can be easily computed using the definition 



of /i in terms of h (from equations ( [3.271 ) and ( |3.28|) ) 



^<^Om^ (413) 

(u = k/a, oj q = q/a, uj c = H c a c /a is the maximal amplified frequency and c^o = 3.1 x 
10~ 18 h wo Hz is the present value of the Hubble parameter; we used that, in our case, 



a(r] c ) ~ b(r] c ) ~ 1). Keeping only the leading terms for kr] c < 1 in eq. (|4.11| ) and (|4.12 ) the 
power spectrum ( 4.13Q is: 



it-,/ v 1 _i , w , 1 , H c vi / ,1, . . ,uj a .,_i 

Sh(uj) ~ z e q i(_ 2(_ 2 — 2(smh7r(^ s, a; > w, 

~ ^^(-)(^)(^)^(sinhvr(^))-i, uj<u q (4.14) 
dJ c Mp Mp lo c 



(where we used that T(iq7] c )T(—iqr] c ) = Tr/(qr] c sinh qr] c ) |p^| ). Since we assumed that the ra- 



diation starts dominating suddenly after — i] c , we can estimate that uj c = 10 11 y H c /Mp Hz, 
assuming that the evolution is adiabatic ( if the evolution is not adiabatic and entropy is 
produced at some stage this result could be slightly modified but for our illustrative purposes 



it is not crucial |37j [ see however for a more quantitative analysis, in a more specific four 
dimensional model]). In order to compare the power spectrum with the phenomenological 
bounds which could constrain the parameter space of our naive model it is useful to work 
with the fraction of critical density stored in the gravity wave background per logarithmic 
interval of uj: 

1 dp GW _! fu q \ ( H c \ 3 2 / u N - 



(4.15) 

Using again the explicit expression for the Bogoliubov coefficients ( |4.11| ) -( [4.12 ) in the two 
different regimes we get: 

&WM - z dec(—) (tt) (sinhvr(^))- 1 , u > u q 



\UJ C J \MpJ ' 'LU C 

(H c /Mp measures how far from the Planck scale the compactification occurs and oo q /u c 
estimates the typical frequency of the internal excitations u q evaluated at the beginning of 
the radiation epoch (77 = — r) c ) with respect to the maximal amplified frequency uj c ; notice 
also that since uj c is the maximal amplified frequency \uo q /uo c \ < 1). While the amplitude 
of the spectra are characterized by the two dimension- less quantities x = log 10 (u q /u c ) and 
y = log 10 {Hc/Mp), the spectral slope is instead fixed by the background evolution and can 
be also more difficult to estimate in a different model of dimensional decoupling. In our 
case the parameter space can be constrained by the observations and since the spectrum 
is increasing in frequency we would expect that the bounds coming from the large scales 
like the COBE bound g7J] (Q gw (uj) < 7.1 x KT 11 for u < uj < 30u ) and the pulsar 
bound || (n GW < 10" 8 at 00 ~ 10 8 Hz) will be less constraining than the bounds arising 
from nucleosynthesis f|9[ (/ dlnujflGwi^, Vo) < 0-2 Qyivo) ~ 1CT 5 , [^7(77) is the fraction of 
critical energy density present in form of radiation, at a given observation time 77]) or from 
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the critical energy density (Q(cj) < 1 for all the frequencies). In particular from ( |4.16|) we 
can find that the COBE bound is satisfied either if y^0.6 x + 35 — 21og 10 /ii o (provided 
wcobe — 10 uj < uj q ) or if y<0.6 x + 8 — 3 log 10 hioo (for ujcobe > oj q )\ the pulsar bounds 
are satisfied either if y^0.3 x + 35 (for oj q < u p ) or if y<2 x + 72 {uj q > uj p ). The COBE and 
the pulsar bounds are less constraining, while the critical density and the nucleosynthesis 
bound combined together give y< — 0.6 x — 0.6 (for uj < to q ) and y^0.3 x — 0.6 (for uj > uj q ) 
which is compatible with H c < 10 -1 Mp and uj q < uj c . 

A similar analysis can be performed in the case of the model ( |2.9|) . The evolution equation 
(4.2) is then (for 77 < — rj c ): 



11 



l a I „ „ ., , , .q 2 r], 



4 z Az 2 



M = 0, z = 2ik(r]-r] c ), a = (4.17) 



which is formally equivalent to the radial Schroedinger equation for the problem of the 
Coulomb diffusion and which can be easily solved in terms of Confluent Hypergeometric 
functions 

fi(kv,q) = e- ik ^^2k(rj - Vc )U(~ - z^, 1, 2ik(rj - rj c )), 77 < rj c (4.18) 

(U is the Kummer function defined with the conventions of [p2]| ; for 77 — > —00 the solution 
behaves like a positive frequency mode , but does not define, asimptotically, a Bunch-Davies 
adiabatic vacuum). As in the previous example we have to match the solution ( 4.1 8|) valid 
for 77 < — 77 c with the solutions (|4.9|) and (|4.10|) valid for 77 > —i] c . The result of this procedure 
will give the Bogoliubov coefficients describing the mixing positive and negative frequency 
modes: 



iiVcY 1 A ,, x! ,2^ 2,v + 2 v 



1kr) c U{~-i^-, 1, -Aikrj c )), k 2 >q 



cx2 



2 /■•'/,. l/ r 

c = e ^[{^/^-y/2k^)U(^-i^-,l, -Aik Vc ) 

+ f-^Ui 1 - - 1, - Aik Vc )), k' < q 2 (H±^h) 2 (4.19) 

Using the small argument limit of the Kummer functions we can compute the normalized 
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spectral amplitude: 



uj c J \MpJ \oj, 

Ship) = ^(smh^rU^W^y^), c< Wj (4.20) 



u; c \MpJ \MpJ \co, 
and the spectral energy density distribution in critical units: 



/ CJg \ u / iJ c \ 8 ( UJ \ 3 , a,' 



- <S)W (siDh <; )r1 © 4 w<w « (4 - 21) 

For uj > uj q the slopes of the spectral energy density distribution (and of the related spectral 
amplitudes) agree with the result previously obtained for the spectra of gravity waves pro- 
duced during a dilaton driven phase , neglecting the internal Laplace-Beltrami operators. 
The amplitudes are instead different due to the presence of uj q . Our partial conclusion is 
that to neglect the internal Laplacians is a good approximation for the slopes of the spectra 
(for uj > uj q ) but not for their amplitude. For uj < uj q both the slopes and the amplitudes 
of the spectra are affected by the presence of the internal Laplacians which cannot be over- 
looked. Since the spectra are increasing in frequency we will keep the most stringent bound 
which comes from nucleosynthesis and which gives, if applied separately in each of the two 
branches, y< — 0.7 x — 0.2 (for uj > uj q ), and y< — 1.1 x — 0.25 (for uj < uj q ). 

The obtained spectral amplitudes ( |4.14| )-( |4~2"0] ) are quite different since they are produced 



by two different background geometries, but the spectral slopes are exactly equal in spite 
of the differences in the solutions ( |2.5| )-( f2~8D . More specifically we obtained "violet" type of 
spectra (5h ~ u/u c for u > u q and 5h ~ (uj/uj c ) 1 I 2 for uj < uj q ) which are a common feature 
of the contracting backgrounds also in the isotropic case [|39|| . This apparent puzzle is due 



to the fact that a~b% ~ ^Iv/Vcl for ( |2.5| ) and ( |2.8| ) in arbitrary number of internal and 
external dimensions. 

We would like finally to point out that the graviton production due to the transition from 
the radiation dominated stage to the matter dominated stage should also be included. This 
further amplification will modify the low frequency tail of the spectrum (10~ 18 Hz < uj < 
10~ 16 Hz). The qualitative aspects of our analysis show that the presence of the internal 
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gradients introduces in the spectral amplitude a new parameter related to the frequency of 
the internal excitation which can be interestingly constrained by the bounds usually analyzed 
in the context of the stochastic gravity-wave backgrounds. 



5 Growing solutions for the scalar modes 



The tensor inhomogeneities can be treated perturbatively keeping track of the internal Lapla- 
cians also because they evolve logarithmically outside the horizon. The situation changes in 
the case of the scalar inhomogeneities because, as we explicitly pointed out at the end of Sec. 
2, the growing solution increases much faster than a logarithm and a perturbative treatment 
become inappropriate after some time. The hope would be in this case that the growing 
mode appearing in the multidimensional case could be consistently gauged down as happens 
for the dilaton driven solutions in the (3 + l)-dimensional case|2^]. In this section we do not 
want to address specifically the problem of the growing modes in an anisotropic manifold 
but we want to show how the problem can be consistently formulated in the presence of 
the internal Laplacians and for this purpose we will study the 10-dimensional dilaton-driven 
solutions ( |2.8|) . Using the all set of equations ( |3.13|) - (|3~T6| ) and equations ( |3.21| ), (|3.22|) , 
( 3.23Q the Fourier modes of the longitudinal fluctuations A can be expressed in terms of 
the Fourier modes of V and W 

11/2 



(k 2 + —q 2 )V(k,q,7 1 ) 



n(n + d- \)HTp' 

3l 2 D <p'H 

[(d-l)H + nF] { a^tf. 

6l 2 D {d- 



Ql 2 D (d-l) 
n(n + d — I) 
V(k,q,r])Y 



W(k,q,T}) 



rap 



lV2 



n{n + d 



(k 2 + -q 2 )A(k J q } 7 1 ) 



(d - l)H + nT 

nT [(n + d- IV] 
d-l\(d- \)H + nT\ 

3l 2 D <p' ( V V 



'W{k,q,r)Y 



Zl 



(5.1) 



6l 2 D (d-l) 
n(n + d — 1) 



nl/2 



'W(k,q,vY 



Zl 



(5.2) 



(d- 1) Va^fet/ 

where q, k have to be considered both scalar eigenvalues (while in the previous section k 
was labeling the eigenvalues of the tensor Helmotz equation). If k 2 > q 2 a 2 /b 2 the evolution 
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equation for V(k,q,rj) and W(k, q, rf) will be (from fl3.21|) , (|3. 22 



) through (|3723|) ) 



V" + 



k 2 + 



4(r/ c - r]f 



V = 0, W + 



k 2 + 



4(r/ c - 77) 5 



W = 



(5.3) 



whose solution is exactly identical to (4.7), provided we choose, for rj — > — oo, the Bunch- 
Davies vacuum as initial condition. In the limit A; 2 < q 2 a 2 /b 2 eq. ( p.21|) -( |3.22|) and (|3.23| ) 
will give instead 



Vcq I 



(rjc-r)) 4(r/ c -r/) 2 
with solution |52| 



V = 0, W" 



V(fc,g,r/) = y/ric- , n H o ) {QVVc-v), W(k,1,v) = VVc ~ V h q Z> (<lVVc ~ v) 



(r/ c - 7]) 4(r) c - rjf 



W = (5.4) 



(5.5) 



The typical amplitude of the longitudinal fluctuations over a length scale k 1 is \5^(k, q,rj) \ = 

d—1 n d-\-rL — 1 

l D k~q2^(k,q,ri) (with Z D = M p 2 ). In the particular case of the 10-dimensional model 
( p.9|) we have \S^(kq, r/)| = l D k 3 ^ 2 q 3 ^k, q (k, q,T]) and from (|5.1|) with the use of ( [B~4l) we 
obtain for A; 2 > q 2 (a/b) 2 



\6*\(k,q,7j) 



He 



\3/2 



(fcr/) 2 



(5.6) 



(we used that Z\ = J\^-\ and Hq (z) ~ \nz). Since A; c ~ l/r/ c we find that \5^(k,q,rj)\ < 1 
on scales k~ x such that \r]/r] c \ > (H c / M P ) 2 (u q / u c ) 3/2 (k c / k) 1/4 . From eq. flOl ) using (|5~5| ) 
we obtain, for k 2 < q 2 (a/b) 2 



15*1 



'A) 3/2 1'k 

k r I 



(6.7) 



,a; c y \k c j \r] / 

which implies that the perturbative approach is only reliable for conformal times 
\ti/ Vc \ > (H c /M P )\co g /co c )(k/k c ) 3 / 2 . 

The presence of the internal gradients slightly changes the quantitative estimates but 
does not change the nature of the growing mode problem. If this is the situation the scalar 
fluctuations will become very soon critical and will then enter in a true non-perturbative 
regime. This apparent contradiction among the behaviour of the tensor inhomogeneities and 
the behaviour of the scalar ones might be the connected with our perturbative technique. 
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The fluctuations which we are discussing now are gauge-invariant only for infinitesimal co- 
ordinate transformations while some quantities which are invariant to all orders would be 
more suitable for the study of fluctuations which are growing outside the horizon. It can 



be actually shown [24] that in the (3 + 1) isotropic case the linearized variables describing 
the scalar and tensor fluctuations in a fully covariant and gauge-invariant approach |fLl| 
grow only logarithmically outside the horizon. Unfortunately a fully covariant and gauge- 
invariant formalism is only formulated in the case of a homogeneous and isotropic manifold 
and it seems to be quite complicated to formulate it for an anisotropic manifold. Within the 
linearized theory discussed in this paper it would be anyway interesting to understand if it 
is possible to "gauge-down" the growing mode solutions arising in a anisotropic background 
perhaps with a suitable generalization of the "off-diagonal" gauge |5(| to the case of an 
anisotropic background. 



6 Conclusions 

We discussed the treatment of the scalar and tensor inhomogeneities in an anisotropic back- 
ground manifold undergoing dimensional decoupling. 

We showed that it is possible to study the evolution equations of the scalar modes by 
completely fixing the coordinate system with a suitable gauge choice which reduces, in 
the isotropic case, to the well known conformally newtonian gauge often employed in the 
analysis of the density fluctuations in the context of the inflationary models driven by a scalar 
field or by perfect fluid matter. The coupled system of second order differential equations 
describing the fluctuations of the scalar field and of a homogeneous and anisotropic manifold 
was diagonalized in terms of two scalar variables which become, in the absence of internal 
dimensions, the normal modes of oscillation obtained by perturbing scalar-tensor action to 
second order in the amplitude of the fluctuations. The evolution equations for the amplitude 
of the scalar perturbations (depending on the internal and external coordinates) were also 
explicitly solved in a particular 10-dimensional background motivated by String Cosmology. 
The scalar spectral amplitudes grow outside the horizon faster than the tensor amplitudes 
and the dependence on the internal coordinates does not change drastically this situation 
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unless the present value of the typical frequency of the internal dimensions would be much 
smaller (by several orders of magnitude) than the maximal amplified frequency. 

The scalar and tensor modes are coupled because the tensors defined in the internal 
manifold are scalar eigenstates of the Laplace-Beltrami operators defined in the external 
manifold. In order to study the connection among the occurrence of a process of dimensional 
reduction and the present energy distribution of cosmic gravitons we discussed explicitly 
two oversimplified toy models of dimensional decoupling. Provided we keep track, in the 
perturbed amplitudes, of the dependence on the internal coordinates the spectral energy 
distribution will be a function of the curvature scale which the compactification occurs and 
of the typical frequency of the internal excitations. We found for both the models a two 
branches violet spectrum. Since the power spectra are increasing in frequency the most 
significant bounds come from small wave-length and constrain the curvature scale to be 
H c < 10 -1 Mp (provided u q < u c and with cu c ~ lO u (H c /M P ) 1/2 Hz). At the same time 
the violet spectra are less constrained at low frequencies by the large scale observations and 
in particular by the CO BE bound. The estimates presented in this paper suggest that for 
uj > uo q the slopes of the spectra (but not their amplitudes) can be reliably computed by 
neglecting the internal Laplacians. On the contrary the presence of the internal Laplacians 
affects decisevely the spectral slopes (and amplitudes) for uj < uj q . Our considerations can 
be also applied to more realistic models of dimensional reduction (with the unavoidable 
help of numerical techniques) in order to discuss the back-reaction problems which could 
eventually lead to the isotropisation of the original background model. It might also be 
of some interest to deepen the possible phenomenological signatures of the scenarios of 
dimensional decoupling and their relevance for the formation of a stochastic gravity-waves 
background. We want finally to stress that even though the models analyzed in this paper 
are quite simplified the perturbative techniques which we introduced are more general. The 
open problem which emerges also from our discussion is of course to understand if a viable 
multidimensional cosmological model, free of the well known problems mentioned in the 
introduction, exists at all. String theory seem to be a very good candidate for this purpose 
and it is very tempting to speculate that the same mechanisms leading, in principle, to 
a graceful exit in four dimensions ||25|| could also operate in order to stabilize the internal 
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dimensions producing, ultimately a completely isotropic universe. It could also be possible 
that classical field configurations (like Dirac monopole configurations polarized along the 
internal dimensions) can offer suitable mechanisms for the stabilization of the internal space 
43| and in this directions the work is still in progress [51 . 
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